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Linofrard's  dielectric  theory  of  energy  loss  is  applied  to  the  stopping  of  energetic  ions  in  solids. 
The  Lindhard  formalism,  in  combination  with  the  local  density  approximation  and  the  concept 
of  effective  charge,  is  used  to  calculate  the  stopping  powers  of  heavy  ions  in  many  targets  and 
over  a  wide  range  of  projectile  energies.  Theoretically  calculated  stopping  powers  are  compared 
with  available  experimental  data. 
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INTRODUCTION  9 

For  many  years,  the  stopping  of  energetic  ions  in  matter 
has  been  a  subject  of  great  interest  to  theoretical  and  experi¬ 
mental  physicists.  Bohr,1'1  Bethe,4-6  and  Bloch’  *  were  the 
pioneering  theoretical  investigators  of  charged  particle  stop¬ 
ping  in  matter.  Some  30  years  after  the  development  of 
Bethe-Bloch  theory,  Lindhard’  presented  an  elegant  dielec¬ 
tric  formulation  of  stopping  theory  for  a  charged  particle 
penetrating  a  free-electron  gas. 

The  Lindhard  dielectric  formalism  probably  stands  as 
the  most  complete  “first  principles”  theory  of  energy  loss  to 
date.  However,  the  application  of  Lindhard’s  theory  to  real 
solids  becomes  a  formidable  task  since  the  direct  numerical 
integration  of  the  Lindhard  stopping  equation  is  slow  and 
cumbersome  due  to  the  presence  of  singularities  in  the  inte¬ 
grand.  This  paper  represents  the  first  attempt  to  apply,  with 
precision,  the  Lindhard  theory  of  stopping  to  actual  solids. 

The  original  work  of  Lindhard  was  extended  by  Lind¬ 
hard  and  Winther'0  (L-W)  who  obtained  an  analytic  expan¬ 
sion  for  the  stopping  power  in  the  limit  of  high  and  low 
projectile  velocity  relative  to  the  Fermi  velocity  of  an  elec¬ 
tron  gas;  the  Bethe  stopping  power  (with  a  redefined  mean 
ionization  potential)  plus  shell  corrections  were  obtained  in 
the  high-velocity  limit  whereas  a  linear  dependence  of  veloc¬ 
ity  with  stopping  power  was  shown  at  low  velocity. 

Bonderup"  utilized  the  L-W  expansions  to  include  the 
local  density  approximation  of  Lindhard  and  ScharfT’  (L-S) 
so  as  toextract  the  Bethe-Bloch  terms,  commonly  referred  to 
as  the  mean  ionization  potential  and  the  shell  corrections.  In 
the  local  density  approximation,  each  volume  element  of  the 
solid  is  considered  to  be  an  independent  plasma  of  uniform 
density,  so  that  the  total  electronic  stopping  power  is  ob¬ 
tained  by  spatially  averaging  a  specific  stopping  power  over 
the  target  electron  charge  distributions.  Bonderup,  in  his 
calculations,  uses  the  Lenz-Jensen  atomic  charge  distribu¬ 
tion  to  describe  the  target  electrons.  His  calculation  was  sub¬ 
sequently  refined  by  Rousseau,  Chu,  and  Powers11  (RCP) 
who  utilized  Hartree-Fock-Slater14  (HFS)  charge  distribu¬ 
tions  to  describe  the  spatially  varying  electron  density,  there¬ 
by  predicting  the  presence  of  oscillations  in  the  stopping 
power  when  plotted  versus  target  atomic  number. 


In  this  paper,  we  evaluate  the  Lindhard  dielectric  the¬ 
ory  of  stopping  power  within  the  local  density  approxima¬ 
tion  and  the  concept  of  effective  charge.  A  comparison  of 
Lindhard's  stopping  powers  with  the  experiment  for  0.4- 
MeV/amu  heavy  ions  in  solid  targets  is  given;  comparative 
agreement  between  theory  and  experiment  over  a  wide  range 
of  projectile  energies  is  discussed. 


LINDHARD’S  THEORY  OF  STOPPING  POWER9 


For  an  ion  of  charge  ze  moving  with  velocity  inn  a 
medium  of  uniform  density/),  the  energy  loss  due  to  electron 
excitation  can  conveniently  be  written  in  the  form 
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where  L  is  the  stopping  number  and  m  is  the  mass  of  the 
electron.  In  the  dielectric  formalism,  L  is  written  as 
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and  e(k,a))  is  the  wave  number  and  frequency  dependent 
longitudinal  dielectric  constant. 

Lindhard  obtained  the  dielectric  constant  for  a  free- 
electron  gas  within  first-order  perturbation  theory  as 
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where  E„  and  k„  represent  the  energy  and  wave  vector  for  an 
electron  in  the  nth  state  and  y  is  a  small  positive  value.  For  a 
free-electron  gas  at  zero  temperature, 


F(£„)=l  for  E„<Eh, 
F  (E„ )  =  0  for  E„  >  £,,, 


(5) 


•Present  address:  Physics  Dept,,  Univ.  of  Chicago,  Chicago,  Illinois  60637  and  the  summation  in  Eq.  (4)  can  be  evaluated  to  get  .  , 
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FIG.  I .  Plot  of  variation  in  stopping  number  [Eq,  {II))  with  electron  density 
for  select  projectile  energies. 
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FIG.  2,  Comparison  of  HFS  radial  charge  density  calculated  by  F.  Herman 
and  S.  Skillman  (Ref.  14)  and  solid-state  spherically  averaged  radial  charge 
densities  by  Moruzzi,  lanak,  and  Williams  (Ref.  16). 
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Here  the  variables  k  and  co  have  been  replaced  by  reduced 
variables 

z=  (k/2k,  )  and  u  =  {a)/kvh  ),  (8) 

where  k,  and  v,  are  defined  by 

Eh  =  Jmtij  =  (/»  2k \/2m)  =  (A  :/2m)(3irp)2/\  (9) 
In  addition,  the  parameter  j-2  in  Eq.  (6)  is  defined  by 

X2  =  (Vu/rrvh),  (10) 

where  v„  =  (e2//t ),  the  Bohr  velocity.  When  the  dielectric 
constant  of  Eq.  (6)  is  substituted  into  Eq.  (2),  one  obtains 
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This  expression  for  L  in  Eq.  (1 1)  is  the  Lindhard  stop¬ 
ping  number.  In  this  paper,  we  make  use  of  stopping  num¬ 
bers  obtained  from  a  numerical  integration  of  Eq.  (11)  as  a 
continuous  function  of  electron  density  (or  y  ’)  over  a  wide 
range  of  projectile  velocities  (or  energies).  While  the  detailed 
results  of  the  calculation  can  be  found  elsewhere,1'  some 
typical  illustrations  of  the  variation  in  stopping  number  with 
electron  density  and  energy  can  be  found  in  Fig.  1. 

985  J  Appl  Phys  .  Vol  51.  No.  2.  February  1980 


LOCAL  DENSITY  APPROXIMATION 

In  the  local  density  approximation,12  each  volume  ele¬ 
ment  of  the  solid  is  considered  to  be  an  independent  plasma 
of  uniform  density  p,  so  that  the  total  electronic  stopping 
power  is  given  by 

(-  f  )-  ZtfYL'pMLW**  <"> 

where  p(r)  is  the  spherically  averaged  charge  density  of  the 
target  atom.  In  this  paper  we  use  spherically  averaged  solid- 
state  charge  densities"'  in  evaluating  Eq.  ( 1 2).  Figure  2  shows 
a  comparison  of  the  HFS  radial  charge  density  calculated  by 
F.  Herman  and  S.  Skillman"  and  solid-state  spherically 
averaged  radial  charge  densities  calculated  by  Moruzzi, 
Janak,  and  Williams'6;  there  is  a  pronounced  difference  in 
the  spatial  variation  of  the  solid-state  and  atomic  charge  den¬ 
sities  away  from  the  interior  of  the  target  atom  due  to  solid- 
state  bonding  effects.  This  spatial  variation  in  charge  density 
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FIG  4  Comparison  of  spatial  variation  in  the  stopping  number  integrand 
of  Eq  (12),  the  stopping  number  of  Eq.  (12).  and  the  solid-state  radial 
charge  density  for  aluminum  (Ref  16)  with  ion  energy  of  10  000  keV/amu. 


gives  rise  to  a  marked  change  in  density-averaged  stopping 
number.  Inspection  of  Figs.  3  and  4  shows  that  the  stopping 
number  integrand  of  Eq.  (12)  approximately  follows  the  ra¬ 
dial  charge  density  in  spatial  variation.  It  should  also  be  not¬ 
ed  that  the  low-energy  (100  keV/amu)  projectile  stopping 
number  is  influenced  mainly  by  the  outer  shell  electron 
charge  distribution,  whereas  the  high-energy  projectile 
(10  000  keV/amu)  stopping  number  is  dependent  upon  the 
inner  electron  shells  as  well.  This  analysis  clearly  shows  the 
importance  of  using  solid-state  charge  densities  in  the  calcu¬ 
lation  of  stopping  powers  of  solids. 

EFFECTIVE  CHARGE  OF  ENERGETIC  IONS  IN  SOLIDS 

In  the  preceding  discussions,  the  projectile  ion  was  as¬ 
sumed  to  be  a  bare  nuclear  point  charge.  However,  when  the 
projectile  is  an  ionized  atom  with  an  intrinsic  electronic 
charge  distribution,  the  question  of  the  net  charge  of  the 
projectile  arises.  From  the  beginning  of  stopping  power  the¬ 
ory  the  charge  of  heavy  ions  in  a  solid  has  been  poorly  under¬ 
stood.  Bohr'  suggested  that  the  ion  might  be  assumed  to  be 
stripped  of  all  electrons  with  orbital  velocities  less  than  the 
ion  velocity  This  concept  is  difficult  to  calculate  in  detail  for 
heavy  ions  because  of  the  difficulties  in  assigning  velocities 
to  the  highly  distorted  electron  orbitals  of  a  partially 
stripped  ion  in  a  dense  electron  plasma. 

The  active  interest  in  fission  fragments  in  1939-1941 
brought  about  various  attempts  to  predict  heavy  ion  stop¬ 
ping  powers  using  the  concept  of  effective  charge. 1  ’ 20  Within 
this  concept,  the  ion  charge  state  was  assumed  to  be  indepen¬ 
dent  of  the  ion-solid  stopping  interaction  so  that  the  ion- 
target  stopping  interaction  chould  be  represented  by  an  ex¬ 
perimental  proton  stopping  power,  and  the  ratio  of  a  heavy 
ion  stopping  power  to  the  proton  value  (for  the  same  velocity 
in  the  material)  was  identified  as  the  heavy  ion  effective 
charge  (it  was  assuried  that  the  proton  charge  was  unity) 
This  quantity  included  not  only  the  charge  state  of  the  ion. 
hut  also  all  reduction  in  slopping  due  to  screening  effects  and 
other  perturbations  (See  review  by  Brandt  for  theoretical 


details.21)  Remarkably,  this  quantity,  for  all  heavy  ions  in  all 
targets  over  several  decades  of  energy,  can  be  represented  by 
a  simple  empirical  formula  established  by  Northcliffe22  as 

z*/z  =  1  —  exp  (v/v0z2/3),  (13) 

where  z*  is  the  effective  charge  and  v0  =  e2h,  is  the  Bohr 
velocity  (2. 1 88-10®  cm/sec).  This  formula  has  been  reshaped 
by  many  authors  as  thousands  of  heavy  ion  stopping  powers 
have  been  measured.  We  use  the  form23 

z*/z=  1  -  [exp(  —  A  )] 

X  [1.034  -0.1777  exp(  -0.08II4z  )],  (14) 

where 

A  =  B  +0.0378  sin \vB, 

(15) 

B  =  0M6(A0E /m)ul/z2/\ 

where  the  ion  energy  E  is  in  MeV,  and  the  ion  mass  is  in  amu. 
The  expression  is  considered  to  be  accurate  to  about  5%  for 
ions  with  atomic  numbers  6-92  in  solids  with  atomic  num¬ 
bers  4-79,  for  energies  above  0.2  MeV/amu.  This  formula  is 
independent  of  target  and  depends  only  on  the  ion  atomic 
number  and  velocity. 

COMPARISON  OF  THEORY  WITH  EXPERIMENT: 
CONCLUSIONS 

Within  the  local  density  approximation  and  the  concept 
of  effective  charge,  we  have  applied  the  Lindhard  stopping 
theory  for  a  charged  particle  in  a  free-electron  gas  to  the 
stopping  of  heavy  ions  in  solids.  Stopping  powers  [Eq.  (12) 
with  z  replaced  by  z*]  were  calculated  for  a  wide  variety  of 
projectile  energies  ranging  from  100  to  10  000  keV/amu  and 
for  many  targets.  For  projectile  energies  of  4000  keV/amu 
and  above,  the  maximum  error  incurred  in  comparing  the 
calculated  stopping  powers  with  experiment2*  for  many  tar¬ 
gets  is  about  4%.  For  projectile  energies  below  4000 
keV/amu,  the  maximum  error  increases  monotonically  to  a 
value  of  about  30%  for  a  projectile  energy  of  100  keV/amu. 
A  typical  comparison  of  theoretical  and  experimental  stop¬ 
ping  powers  for  many  targets  is  illustrated  in  Fig.  5.  Here  the 
ratio  of  experimental  to  theoretical  stopping  is  plotted  versus 
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target  atomic  number  for  a  projectile  energy  of  400 
keV/amu. 

The  objective  of  this  paper  has  been  to  apply  the  Lind- 
hard  theory  of  stopping  to  real  solids.  Within  the  assumed 
approximations,  the  results  show  that  the  theory  significant¬ 
ly  deviate  from  experiment  in  the  low  projectile  energy  (100 
keV/amu<£<4000  keV/amu)  region.  The  reasons  for  this 
disagreement  may  be  manyfold:  the  use  of  a  free-electrdn  gas 
dielectric  response  function,  the  extrapolation  of  the  uni¬ 
form  electron  gas  model  to  solid  state  systems  via  the  local 
density  approximation,  and  the  use  of  the  effective  charge 
concept  to  reduce  experimental  stopping  powers  to  stopping 
numbers.  In  a  future  paper,  we  plan  to  explore  an  extension 
of  Lindhard’s  theory  to  self-consistently  include  heavy  atom 
projectiles  which  may  abrogate  some  of  the  above 
weaknesses. 
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